CrRb: a molecule with large magnetic and electric dipole moments 
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We report calculations of Born-Oppenheimer potential energy curves of the chromium-rubidium 
heteronuclear molecule (^^Cr*^Rb), and the long-range dispersion coefficient for the interaction 
between ground state Cr and Rb atoms. Our calculated van der Waals coefficient (Ce = 1770 a.u.) 
has an expected error of 3%. The ground state ''E^ molecule at its equilibrium separation has a 
permanent electric dipole moment of de{Re = 3.34 A) = 2.90 D. We investigate the hyperfine and 
dipolar collisions between trapped Cr and Rb atoms, finding elastic to inelastic cross section ratio 
of 10^-10^ 

PACS numbers: 34.20.-b,34.50.Cx 



I. INTRODUCTION 

The creation and manipulation of ultracold polar 
molecules in their ground state is a growing endeavor 
in atomic and molecular physics. Rubidium is the 
workhorse of ultracold atomic physics. It was the first 
atom to be Bose condensed and it remains a favorite atom 
for simulations of correlated and many-body physics, and 
for forming cold molecules. The chromium atom has also 
been Bose condensed and it is the only atomic species 
for which strong dipolar interaction has been observed 
[i|. Dipole moments, by virtue of their long-range 1/R^ 
interaction, where R is the separation between atoms, 
can be tuned by external fields with increasing precision 
to control atomic interactions with applications to pre- 
cision measurements, molecular physics and information 
processing. 

CrRb is a heteronuclear molecule, and may therefore 
possess a sizable permanent electric dipole moment. In 
its ground electronic state, it has magnetic moment of 
five Bohr magnetons (S/is), so it can be magnetically 
tuned. In its most abundant form (84%), ^^Cr has no 
nuclear spin, i — 0, but its fcrmionic isotope ^■^Cr (9.5% 
abundance) has nuclear spin, i = 3/2, that couples to 
electronic spin, s = 3, to produce a number of hyperfine 
levels. Therefore, bosonic-fermionic mixtures of CrRb 
can be formed with large electric and magnetic dipole 
moments, with potentially interesting applications for de- 
generate dipolar Fermi gases and spinor physics. 

There is no spectroscopic information available for 
the CrRb molecule. A two-species magneto-optical trap 
(MOT) for Cr and Rb was reahzed in 2004 g, in which 
some 4 x 10^ ^^Cr atoms and 3 x 10^ ^^Rb atoms were 
loaded. Cr is known to have large inelastic two-body 
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spin-flip losses, so it cannot be maintained in a MOT 
[si y|. Owing to its large magnetic moment, magnetic 
trapping (MT) of Cr is possible Q, which allows the 
atoms to be trapped in its lowest high-field seeking 
state with total electronic spin s = 3 and its projec- 
tion nis = -3. In the Rb-MOT + Cr-MT configuration, 
Hcnsler et al. 2] measured the two-body loss rate constant 
to be /3RbCr ~ 1-4 X IQ-" cm'^/s and /3crRb '- 10"^° 
cm^'^/s, where the former refers to the loss due to the in- 
troduction of Cr into the Rb-loaded MOT, and the latter 
refers to loading of the Cr-MT first. The latter loss coef- 
ficient is about an order of magnitude larger than the for- 
mer loss coefficient, because in the more shallow Cr-MT, 
interspecies dipolar interactions lead more quickly to spin 
depolarizing collisions, hence depletion of the trap. 

In this ab initio work, we analyze the electrostatic in- 
teraction of Rb and Cr atoms in their ground states, 
obtain Born-Oppenheimer (BO) potential ground state 
energy curves, the long-range dispersion Cg coefficient, 
and present a value for the electric dipole moment of the 
CrRb molecule. We use the resulting potential energy 
curves to calculate the elastic and relaxation cross sec- 
tions, hyperfine and dipolar Fano-Feshbach resonances 
which can be used to tune the interactions. 



II. NUMERICAL CALCULATIONS AND 
COLLISIONAL RESULTS 



The calculations of the BO potential energy curves 
and the permanent electric dipole moments were per- 
formed with the multiconfiguration complete active space 
self-consistent field (CASSCF/CASPT2) method [1 §. 
The basis set was of VQZP quality with the primi- 
tives obtained from the relativistic ANO-RCC basis set 
{7s6p4:d3f2glh for Cr and 8s7pAd2flg for Rb)[7i, Q. 
Scalar relativistic effects are included in the calculations 
using the Douglas-Kroll-Hess Hamiltonian, as is standard 
in the MOLCAS software. Two active spaces were used. 
The first comprised the Cr 3d and 4s orbitals and the Rb 
5s, thus seven orbitals with seven active electrons. In the 
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FIG. 1: The '^E+ and ®E+ potential energy curves for CrRb, 
correlating to Cr('^S'3) and Rb(^5'i/2) ground states. 



second set of calculations, a second set of 3d orbitals was 
added to describe the 3d double shell effect: 12 orbitals 
with seven electrons. All calculations were performed 
with the MOLCAS-7 quantum chemistry software ^]. 

There are two electronic states which correlate to Rb 
(6s ^S) and CT{3dHs ^S") ground states, ^E+ and ^T.+ . 
At large separation, R, these curves are well described 
by —Cq/R^, where Cg is the van der Waals (vdW) co- 
efficient. We calculate Cg using the Casimir-Polder [lO'| 
integral 
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where the dynamic polarizability of each atom, a(ia;), is 
given by 
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The calculation of the dynamic polarizability requires 
knowledge of the oscillator strengths /qa between the 
atomic ground state of energy cq and atomic excited 
states of energy e\. The summation is understood to 
also include integration over continuum. Our calculated 
value of Cg for the Cr + Rb system, 1770 a.u., which 
we believe to be accurate within 3%, is obtained from 
the highly-accurate values of the Rb dynamic polarizabil- 
ity at imaginary frequencies |11| and the recently accu- 
rate values for the Cr dipole polarizability [12, [l3| . This 
should be compared to the values of C6(Rb2) = 4691 a.u. 
[13, [Hi, and CeCCra) = 770 a.u., 4]. Both ^2+ and ^1]+ 
curves share the same Cg coefficient. 

In Fig. [TJ we plot the BO potential energy curves for 
CrRb in the ground electronic states. The equilibrium 
separation for the molecule in the ground state ^S^ 
is Re = 6.31 a.u. The value for the dipole moment 
at the equilibrium distance for the ^E"^ ground state 
molecule is de = 2.90 D (1.14 a.u.), which is reasonably 
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FIG. 2: (Color online) Top panel: Elastic cross section for the 
^S"*" and ®S+ curves without shift in the inner wall {s — 0). 
Bottom panel: Effect of the positive shifts on the elastic cross 
section along ^S"*" curve. Potential energy points for which 
R < Re are shifted according to Rs = R+s{R— R^) / {Rt — Re) , 
where s is the shift of the zero-energy classical inner turning 
point, Rt, and Re is the equilibrium distance. Several different 
values for s were used. 



large. The calculated static dipole polarizability for the 
molecule is Q;crRb(0) = 701 a.u. We should mention that 
aRb(O) = 319 a.u. [H] and acr(O) = 86 a.u. [4]. The 
molecular value is 40% larger than the sum of atomic 
polarizabilities. 

The dearth of spectroscopic and coUisional information 
about the CrRb molecule limits what can be extracted 
from our calculations at ultracold temperatures. Never- 
theless, we can inform the discussion by modifying the 
interaction potentials and evidence the changes to the 
phase shifts and/or cross sections at very low energies. 
Our vdW coefficient is sufficiently accurate to describe 
the potentials at large separations. In the short range re- 
pulsive wall region, we shift the potential energy curves, 
according to the prescription in our earlier work jl2| . 
and the results for the elastic cross sections are shown 
in Fig. m Several points can be made: while the cross 
sections are insensitive to the shifts for energies above 
E > 10~^ a.u., they tend to decrease (increase) with in- 
creasing positive (negative) shifts at very low energies. 
However, when the shifts are large enough that an addi- 
tional bound state (see for instance Fig. [3|) is pulled in 
from the continuum, the above trend does not hold. 

III. ZEEMAN CASCADE 

The scattering between Cr and Rb atoms in the pres- 
ence of a homogeneous magnetic field, B, is governed by 
the following Hamiltonian 
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FIG. 3: The dependence of the energy of the last bound state 
in the ^S"*" and *E^ potentials on the shift parameter (s), see 
Fig. [2] for the shift definition. The "resonances" occur when 
an additional bound state is pulled in from the continuum as 
s decreases. 



where the constants /i, fiB and Ahf represent the re- 
duced mass of the CrRb molecule, Bohr magneton, and 
Rb isotropic hyperfine interaction constant of the ground 
state, respectively. The first two terms in the Hamilto- 
nian represent the nuclear kinetic energy of the molecule, 
where I is the rotational angular momentum of the nu- 
clei. The explicit form of Vcs, the operator of the electro- 
static interaction, and Vdip, the operator of the magnetic 
dipolar interaction will be given below. The spin opera- 
tors saIia) and SB(iB) represent the electronic (nuclear) 
spins of atom A (Cr) and atom B (Rb), respectively. For 
^^Cr + ^^Rb system, the electronic and nuclear spins are 
SA = 3, ZA = 0, SB = 1/2 and zb = 3/2 The linear 
atomic Zeeman terms and the isotropic hyperfine inter- 
action for Rb are also given. The i?-independent terms in 
the Hamiltonian ([3]), determine the scattering channels, 
which can be related through a unitary transformation 
1/^) = J2a Uj3a\ct), to the product states of the total elec- 
tronic spin, S' = sa + sb, and nuclear spin, I = XA + iB, 
\a) — \SMsIMi). In this work, we neglect the interac- 
tion of the Rb nuclear magnetic moment with B. 

The total wave function expansion by product states 
of TsMs 1 the eigenfunctions of S^ and Sz , and tjmi , the 
eigenfunctions of /^ and Iz, is: 



^= 2_^ 2_^ GsMsIMiiR)TSMs'''IMi- 



(4) 
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We introduce channel wave functions, GsmsIMi (-R)? that 
describe the motion of the nuclei. Substitution of the ex- 
pansion (HI) in the Schrodinger equation leads to a system 
of coupled differential equations. A further simplification 
is achieved by introducing a set of channel wave functions 
that depend only on R, 
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where R determines the direction of R. The operator V^s 
is defined in terms of the Born-Oppenheimer interaction 
potentials Vs{R) for the CrRb molecule, 



Ves = J2J2 VsiR)\SMs){SMs\ 

S Ms 



(6) 



This decomposition makes the operator V^cs diagonal in 
the basis (gj. 

The magnetic dipolar operator [l5| , in the second rank 
irreducible tensor representation, can be written as 
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where a is the fine-structure constant. The matrix of 
the dipolar interaction ([7]) can be evaluated analytically 
as demonstrated, for example, in [ial3- I^ i^ ^^e only 
term in the Hamiltonian that can couple channel wave 
functions with different rotational numbers, I, according 
to 



A/ = 0, ±2; while 0^0 is forbidden. 



(8) 



A similar rule holds for the change of the total electronic 
spin of two atoms 



AS' = 0,±1,±2. 



(9) 



The dipolar interaction preserves the angular projections. 
Ml, and Ms + mi independently. 

The Zeeman interaction term is 2iiqBMs, where Ms 
is the projection of S on B. The hyperfine interac- 
tion term can be directly evaluated in the coupled basis 
IfAfsFMr), where /a = ia -H sa, /b = ^b -H sb, and 
F = fA + fs, and then through the chain of transforma- 
tions 



IfAfBFMp) ^ \SIFMf) -^ \SMsIMi), 



(10) 



expressed in the \SMsIMjlmi) basis. The hyperfine in- 
teraction matrix elements are calculated as in [ij, [I^ . 

The highest low-field seeking state is chosen as the 
initial channel for the Zeeman relaxation of maximally 
stretched Cr and Rb atoms. In this state, the spin num- 
bers for ^^Cr are / = 3 and nif = 3, while for ^^Rb the 
relevant spin numbers are f — 2 and ruf — 2. The results 
do not change much when we exclude states which rep- 
resent the decrease of the total spin projection, Mp, by 
more than 2 units. Figure 4(a) contains the cross sections 



as a function of energy at three values of B. The cross 
sections display characteristic dipolar shape resonances, 
and roughly correspond to cross sections for Cr-Cr scat- 
tering scaled down 10 times. The inelastic cross section 
contribution from the dipolar interaction scales as [l9[ 
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FIG. 4: (Color online) (a) Energy dependence of inelastic 
cross section for the relaxation of maximally stretched ^ Cr, 
|/=3, r7i/=3), and *^Rb, |/=2,my=2), atoms in the presence 
of magnetic field, (b) Rate constants for elastic and inelastic 
scattering processes presented in (a). 



The Zccman relaxation rates, Fig. |4(b)[ provide a valu- 
able tool in assessing the efficiency of CrRb evaporative 
cooling in a magnetic or optical trap. Ratios of elastic to 
inelastic collisions in the fiK regime of about 100—1000 
times are possible. 

In CrRb, we have an interesting system where 
strong hyperfine and dipole-dipole interactions 
compete to produce rich spectra. When only the 
hyperfine interaction in ^^Rb is considered, the 
scattering length dependence with the magnetic 
field. Fig. [5] (bottom panel) shows a number of 
hyperfine-coupled Fano-Feshbach resonances, where 
the initial state is the lowest high-field seeking state, 
|52Cr : /=3,m/=-3; «^Rb : /=l,m/=l; /=0,m,=0). 
When the Cr magnetic dipole interaction with the spin 
of the Rb atom is "turned on" , the much richer spectrum 
(top panel in Fig. [5]) is obtained, where the additional 
resonances, due to the dipolar interaction, have (s, d) 



wave character. 

In summary, we calculate the ^E+ and ^I]+ BO poten- 
tial energy curves of CrRb in the ground electronic state, 
the long-range Ce vdW coefficient, the static molecu- 
lar polarizability, and the equilibrium permanent elec- 
tric dipole moment. We obtain the spin relaxation rate 
constants for collision of Cr and Rb in a magnetic field 
and investigate the magnetic tunability of the scatter- 
ing length in the presence of hyperfine and dipolar in- 
teractions. With its relatively large electric and mag- 
netic dipole moments, CrRb is a potentially interesting 
molecule for coUisional studies of dipolar molecules in 
combined electric and magnetic fields. With fields at dif- 
ferent angles, it may be possible to tune the two dipoles 
into cooperative or competitive arrangement. 
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FIG. 5: Competing hyperfine and dipolar interactions in the 
I = entrance channel state: |/=3, m/=— 3) for ^^Cr, and 
|/=l,r7i/=l) for *^Rb. Bottom panel: Resonances caused by 
Rb hyperfine interaction. Top panel: Effect of dipolar inter- 
action coupling entrance s-wave channel to the closed d-wave 
channels. 
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